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Abstract 



Based on a perturbative theory of quantum chromo dynamics and 
non-relativistic quark model, associated J/ ip plus open charm photopro- 
C^j '■ duction on charm quarks in a proton via partonic subprocess 7 c — ► J / ipc 

is discussed. It is shown that the value and energy dependence of 
the cross section for such process remarkably depends on the choice 
I of charm distribution function in a proton. In the region of small 

z = Ej/Ej < 0.2 the contribution of the 7c — ► J /ipc subprocess in the 
inelastic J ftp photoproduction spectra is larger than the contribution 
I of the photon-gluon fusion subprocess. At the energy range of HERA 

■ collider intrinsic charm contribution in the total inclusive J/tp photopro- 

duction cross section may be equal to 4% of the dominant contribution 
of photon-gluon fusion mechanism. 

■ 

^ ■ 1 Introduction 

a.: 

^ ■ The study of the charm quark distribution function in a proton taiks a remarkable interest 

I from the point of view of a investigation of the non-perturbative proton wave function Q 

1 as well as from the point of view of a calculation of the intrinsic charm quark contribution 

to the processes of charm particle and quarkonium production in photon-hadron, hadron- 
•h , hadron and hadron- nucleus interactions at high energies |^|, ||. The existing data on 

charm production processes gives us the idea [Q, [I| of a small, but finite (0.3%-0.5%), 
non-perturbative cc-component in proton wave function. It helps to explain some effects 
in charm production processes which are difficult for understanding in the assumption that 
charm quarks were produced only in hard partonic subprocesses Q. 

The J ftp photoproduction process is the best source of information about gluon struc- 
ture function in a proton especially in the region of very small In the deep 
inelastic domain J/tp photoproduction process may be described using so-called colour 
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singlet model [M p|, in which J/ip is produced in 717 — > J/ipg partonic subprocess. Obvi- 
ously, in the proton fragmentation region the contribution of the intrinsic charm in proton 
via 7c — + J I i/jc subprocess is too comparably large Q . The partonic subprocess 7c — > J/ tpc 
gives the main contribution to associated J /if) plus open charm photoproduction and may 
be very clean test for charm quark sea in a proton. 

In Sec. 2 we shall obtain amplitude, differential and total cross section for partonic sub- 
process 7c — > J /tpc using perturbative theory of QCD and non-relativistic quark model. 
The connection of the partonic and measurable cross sections for processes "fp — > J/ipcX 
and 7p — > J/tpX is discussed in Sec. 3. In that part we also present charm quark distri- 
bution functions which are used in calculations. In Sec. 4 we shall calculate z-spectra and 
total cross sections for associated J ftp plus open charm photoproduction process as well 
as the 7c — > J /tpc subprocess contribution to the inclusive J/tp photoproduction cross 
section at high energies. 



2 Partonic subprocess 7c — > J/ipc 

In the lowest order of QCD perturbative theory associated J/tfj plus open charm pho- 
toproduction in 7p-interactions corresponds to partonic subprocess 7c — > J / ipc* It is de- 
scribed by Feynman diagrams which are shown in Fig.l. The quarkonium is represented as 
non-relativistic quark-untiquark bound system in singlet colour state with specified mass 
M = 2m (m is c-quark mass) and spin-parity J p = l~ . 

The amplitude of subprocess 7c — > J/ipc can be expressed in the form: 

M = M 1 + M 2 +M 3 +M i , (1) 



where 



Ma = e^(^^ ( ;_' fc)2 _ m ^ U iq) , (4) 
M 4 = e^U^h^ j^^ &rP^^^Uiq). (5) 
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In these formula: 

P=-^Aej(p + m), 

A = *(0)/Vto, F c = S kr /Vi, k and r are colour indexes of charm quarks, T a — A a /2, 
e q is electrical charge of -quark in units of e. It is well known that ^(O), which is equal 
to J ftp wave function at zero point, can be extracted in the lowest order of perturbative 
QCD from the leptonic decay width of the J/tp: 

H TTl 

We shall put in our calculation T ee = 5.4 KeV fiJJ. 

If we average and sum over spins and colours of initial and final particles, we obtain 
the expression for square of matrix element: 



l^l 2 ^ £ (7) 
j>i=i 



where 

32ir 2 a 2 T ee m 

B~,r. — ~ 



9a 

2 I _ i/< 2 - _ - / 2 



s = s/m , t — t/m , u — u/m 

s, t, u are usual Mandelstam variables and s + i + u = 6m 2 . The explicit analytical formula 
for functions Kij have the following forms: 

K u = -(2st-2s + Pu-AP-8tu + lAt + 7u-106) 

/(4 * (t 4 - Ap + 6? - 4i + 1)) (8) 

K12 = (s 3 - s 2 t - 6s 2 -st 2 - 2siu + 16st - su 2 + 8su 

-28s + P - 6P - iu 2 + 8tu - 28t + 4u 2 - 126u + 276)/ 

(4(S 2 ? - 2s 2 t + s 2 - 2st 2 + Ast - 2s + P - 2t + 1)) (9) 

K 13 = (s 2 t + s 2 + 2sP + stu - 22st + 7su + 16s + 

2t 2 u - 2AP - life + 217£ - 2u 2 + ilu - 511)/(2 * (s?« - 4S 

i 2 - 2stu + 8st + su - As - Pu + At 2 + 2tu - 8t - u + 4)) (10) 
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if 14 = - (s 3 + s 2 i - 18s 2 + sP + stu - 16s* - su 2 - 3su + 
166S - ll* 2 - htu + 115*+ 13m 2 - 65u - 335)/ 
(2 * [Pu - 4* 3 - 3Pu + 12t 2 + 3tu - 12* - u + 4)) (11) 

if 22 = -(S 2 u - 4.5 2 + 2st - 8su + US- 2t + 7u - 106)/ 

(4 * (s 4 - 4s 3 + 6s 2 - 4s + 1)) (12) 

if 23 = -(s 2 * - lis 2 + sP + stu - 16st - 5su + 115s + 

* 3 - 18? - tu 2 - 3tu + 166* + 13u 2 - 65u - 335)/ 

(4(S 3 u - 4S 3 - 3S 2 u + 12S 2 + 3Su - 12s - u + 4)) (13) 



if 24 = (2s 2 * + 2s 2 i2- 24s 2 + sP + stu - 22s* - llsu + 

217S + * 2 + 7tu + 16*- 2u 2 +41u- 511)/(2* (S 2 tu - 4S 2 t 

-s 2 u + AS 2 - 2stu + 8St + 2su -8§ + tu-4t-u + 4)) (14) 

if 33 = -(2S 2 + 2st + 6SH- A&s + Pu- 10* 2 - 

Atu + 74* + 2u 3 - 22u 2 + 53u - 118)/(S 2 u 2 - 8S 2 u + 

16S 2 - 2su 2 + 16su - 32s + u 2 - 8u + 16) (15) 

if 34 = -(2 * (S 2 + st- lls + P + 2tu - 19* + u 3 - 9u 2 + 28u + 11))/ 
(stu 2 — 8stu + 16s* — su 2 + 8su — 16s 

-tu 2 + 8tu - 16t + u 2 - 8u + 16) (16) 

if 44 = - (s 2 u - 10s 2 + 2s* - 4su + 74s + 2P + 6tu - 

46* + 2m 3 - 22u 2 + 53u - 118)/{Pu 2 - 8Pu + 16* 2 - 2tu 

2 + 16*u - 32* + u 2 - 8u + 16) (17) 

The differential cross section for subprocess 7c — ► J/ipc can be written as follows: 

do 1 , 2 



dt 16tt(s - m 2 ) 2 



|M| 2 (18) 
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The total cross section will be obtained after integration over t in limits: 



t> 



2s 



■ [s — 3m 2 ± y/ (s — 9m 2 ) (J — m 2 )] . 



This procedure can be made analytically, we find that: 
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a(jc — > J/ipc) 



13- 



16tt(s — rn 2 ) 2 



(19) 



j>i=i 



where t r , 



c/rn 2 and t r , 



t-min /m 2 - The explicit expressions for functions 



Hij(s,t) are more unwieldy than for functions Kij and the FORTRAN expression for 
can be obtained by E-mail from the author on request []. 

We shall also calculate the dominant contribution in total and differential inclusive J/tp 
photoproduction cross sections from 717 — * J/ipg partonic subprocess. Here we present 
main formula for this one without discussion (see, for example j7j). 

First the partonic differential cross section is given by: 



where 



Here: s + i + u — Mj, i r 



^( 7 <? -> J/Vfl) = B^ g MjF{s,% 
at 



_ 87TQ 2 r ee 

19 - 3aMj ' 

l/f2\2 



s 2 {s - Mj) 2 +t 2 {t- Mj) 2 + u 2 (u~ MjY 



( S -Mj) 2 (t-Mj) 2 (u-Mj) 2 



— £min — 



Mj 



The total partonic cross section reads: 

(7(75 -> J/i>9) = B ig Mj$(s) 

where 

$(.s) 



M 2 



21n(s/M 2 ) 



(s + Mj) 2 
2(3 + Mj) 



sMj 



Mi 



41n(s/M 2 ) 



s 2 Mj{s-Mj) s{s-Mj) 2 ' 



(20) 



(21) 
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3 Associated J ftp plus charm photoproduc- 
tion on proton 

Let us consider the kinematic for process 7p — > J/ip + c + X in the rest frame of the proton 
(the lab. frame). Variables of partonic subprocess 7c — * J/ipc and variables describing 
measurable process are connected as follows: 

s = xs + to 2 , t = 5m 2 — xzs, u — — xs(l — z), 

where (in lab. frame) s — 2m p E 1 , z — Ej/E 7 , a; is a momentum fraction of charm quark 
in the proton. In the general factorization approach of QCD the measurable cross section 
a and partonic cross section a are connected by the following expressions: 



cr( 7 p^ J/tJjcX) = dxC p {x,Q 2 )a{~fc^ J/tpc), (22) 

where Xmin — 8m 2 /s, Q 2 = Mj and 

J/ipcX) =s dxxC p (x,Q 2 )-^(jc^ J/tfjc), (23) 
dz J dt 

where the region of integration over x is defined by the condition: t m i n < 5m 2 — xzs < t max . 

At present the direct experimental information about charm quark distribution func- 
tion in a proton in wide region of x is practically absent. Existing parameterizations of 
C p (x, Q 2 ) are very different. For comparison we shall use in our calculation "hard" scaling 
parameterization Q: 

C p (x,Q 2 ) = C p {x) = 18a; 2 [(l-x)(l + 102: + a: 2 )/3 + 2a;(l + a;)lna;] (24) 

and "soft", based on perturbative QCD, parameterization at the scale Q 2 = Mj: 

xC p {x, Q 2 ) = (s - s c ) a (l + Bx)(l - x) D cxp (-E + ^ E 1 s b / x)^j , 

where 

o=1.01, 6 = 0.37, s c = 0.888, 5 = 4.24- 0.8045, D = 3.46 + 1.076s, 
E = 4.61 + 1.49s, E' = 2.555+ 1.961s, fi 2 = 0.25 GeV 2 , A = 0.232 GeV : 



(25) 
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/ ln(Q 2 /A 2 ) \ 
W( M 2 /A 2 ) J 



Note that mean value of the proton momentum, which is carried by charm quarks, 
is equal to 0.3% in the case of parameterization [EJ and approximately 0.5% in the case 



of parameterization 11 , which does not contradict data from EMC Collaboration on 
Fi(x,Q 2 ) OJJ. 

These parameterizations for C p {x, Q 2 ) have very different physical interpretation. The 
scaling parameterization 0], so-called "intrinsic", was obtained assuming existence of a 
small, but finite, non-perturbative charm component in proton wave function. It is neces- 
sary for description of open charm and J /if) production total cross section and xf spectra 
at Xf — > 1 in hadron-hadron and hadron-nucleus collisions. On the contrary, the param- 
eterization Jll]|, so-called "extrinsic", strongly depends on choice of scale Q 2 , because it 
was obtained in perturbative QCD approach assuming that at Q 2 < charm quarks 

in the proton are absent and they are generated in QCD cascade only at large Q 2 . 

In Fig. 2 the x dependence of the charm distribution function is shown at Q 2 = M 2 . 
As it will be discussed later different x- dependences of the charm distribution functions 
give us very different predictions for total J / ip photoproduction cross section as a function 
of energy. 



4 Results and discussion 

Result of calculation of total cross section as function of photon energy E 1 for associated 
J/ip plus open charm photoproduction in approach discussed above is shown in Fig. 3. 
"Intrinsic" parameterization Q predicts largest value of cross section at energies E 1 < 60 
GeV, but this cross section rapidly decrease as energy growth at E 1 > 100 GeV. On the 
contrary, the cross section, which was calculated using "extrinsic" parameterization |TT) , 
monotonously increase from 510 -3 nbn at the energy E 1 = 50 GeV to 0.4 nbn at the 
energy E 1 = 5 TeV and at the energy range of ep-collider HERA {\/s — 200 GeV) it is 
approximately equal to 0.8 nbn (Fig. 4). Such a way we have sufficiently large measurable 
cross section of the J/ip plus open charm production for direct experimental investigation 
of the x-dependence of charm quark distribution function in a proton. 

The study of the J/ip photoproduction in jp- interactions concentrated on inelastic 
J ftp production as a tool to obtain information on the gluon distribution function in a 
proton ||. That is why we have calculated the contribution of the 7c — > J/ipc partonic 
subprocess in the total inclusive J ftp photoproduction cross section. Fig. 4 shows result of 
calculation this contribution using "extrinsic" parameterization as well as contribution 
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of the dominant photon-gluon fusion subprocess as functions of \fs . It is known that 
7<7 — + J/ijjg partonic subprocess gives contribution in the total inclusive J/ip photopro- 
duction cross section approximately equal to 50% at the energy \fs lp — 200 GeV Q . Our 
calculation shows that total contribution of the subprocesses 7c — > J/ipc and 7c — > J/ipc 
at the energy y/s = 200 GeV is equal to 4% of the photon-gluon fusion subprocess con- 
tribution or 2% of the total contribution of all mechanisms. The contribution of the charm 
quarks in J/ip photoproduction cross section may be significant and approximately equal 
to contributions of elastic (5%) or diffractive (2%) mechanisms ||. But at smaller ener- 
gies y/s — 10 — 20 GeV charm quark and untiquark contribution is 15% of photon-gluon 
fusion contribution in the case of parameterization 0] for charm quarks in the proton. 

The contribution of the proton intrinsic charm in J/tp photoproduction must be very 
large in the region of proton fragmentation that is at small z = EjjE~ i in the proton 
rest frame. Figs. 5 and 6 show results of calculation for contributions 7c — > J/ipc and 
7<7 — > J/ipg partonic subprocesses in z— spectra of the J/ip photoproduction at energies 
\fs — 14.7 GeV and y/s — 200 GeV, correspondingly. The charm quark contribution 
is larger than photon-gluon fusion subprocess contribution at z < 0.2 for both energies. 
At the relatively small energy (y/s — 14.7 GeV) the parameterization |jj gives more 
large contribution. At the energy \fs lp — 200 GeV the contribution of parameterization 
Q is largest only for very small z < 10~ 3 , in the region 10~ 3 < z < 0.2 the contribution 
of parameterization pd| is dominant. 

In conclusion we note that the comparison of our results with data didn't made because 
data usually have been obtained in fixed kinematical region of variables z and pt0- In 
contrary in our approach total kinematical region (It is bounded by conservation laws) of 
these variables have been took into account. But all calculations were made in the same 
kinematical approximation and we can to compare relative contributions of the different 
mechanisms with data. The calculations for total cross section, z and pt spectra of the J/ ip 
photoproduction via partonic subprocess 7c — > J/ipc in fixed kinematical region (HERA 
collider, for example) will be present in future publications. 

I would like to thank Prof. N.Zotov and Dr. A.Martynenko for useful discussions 
and particularly Prof. A.Likhoded for remarkes. The work was supported in part by the 
Russia Foundation for Basic Research under Grant 93-02-3545. 
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Figure captions 

1. Diagrams used to describe the partonic subprocess 7 c — > J /ipc. 

2. The charm distribution function in a proton versus x at Q 2 — Mj. The curver 1 
corresponds to parameterization fiL the curve 2 - fTl{| . 

3. The cross section of the process 775 — > J/ipX as a function of photon energy E 1 . 
Curves as the same Fig. 2. 
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4. The cross section of the process jp — ► J/ipX as a function of y^s . The curver 1 
corresponds to subprocess 7c — ► J/ipc (parameterization for C p (x,Q 2 ) is taken from 
p| ), the curver 2 - 73 — ► J/"0g (parameterization for G p {x,Q 2 ) is taken from ]Tl[ | 
too). 

5. The z- spectrum of J/0 in jp- interaction at y/s — 14.7 . Curvers 1 and 2 - 
contributions of 7c — > J/ipc subprocess used parameterizations Q and H], corre- 
spondingly. The curver 3 - contribution of photon-gluon fusion mechanism. 

6. As the same Fig. 5 at y/s = 200 GeV. 
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Abstract 



Based on a perturbative theory of quantum chromodynamics and non-relativistic quark 
model, associated J/tp plus open charm photoproduction on charm quarks in a proton 
via partonic subprocess 7c — > J/ipc is discussed. It is shown that the value and energy 
dependence of the cross section for such process remarkably depends on the choice of 
charm distribution function in a proton. In the region of small z = Ej/E 1 < 0.2 the 
contribution of the 7c — > J/ipc subprocess in the inelastic J/ip photoproduction spectra is 
larger than the contribution of the photon-gluon fusion subprocess. At the energy range 
of HERA collider charm quarks contribution in the total inclusive J/ tp photoproduction 
cross section may be equal to 4% of the dominant contribution of photon-gluon fusion 
mechanism. 
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